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I. INTRODUCTION 

There is a considerable interest to understand the gravitational dynamics of asymptotically anti de Sitter (AdS) 
space-times, stimulated to a large extent by the AdS/CFT correspondence, but independently of that, the problem 
is of sufficient inherent interest on its own In the ground breaking work |lfll |. the time-evolution of a free, 

massless scalar field coupled to Einstein’s gravity has been investigated in 3 spatial dimensions in asymptotically 
AdS space-times, with the result that from a large class of smooth initial data black holes form, indicating that AdS 
is unstable. By now a large body of rather convincing numerical and perturbative evidence has accumulated that 
asymptotically AdS s pac e-times (having the same conformal boundary that AdS) are in fact unstable with respect to 
black hole formation [lll - l^ . The fact that arbitrarily small perturbations of an AdS space-time lead generically to 
black hole formation is in sharp contradistinction to the Minkowski or de Sitter cases which are stable. The instability 
uncovered by Ref. manifest itself by the fact that more and more energy gets concentrated in the same spatial 
region by weak turbulence. The existence of such an instability is clearly related to the peculiar causality structure 
of asymptotically AdS space-times, which are not globally hyperbolic, ie. there is no Cauchy hypersurface in them. 
As a consequence, the prescription of initial data on a space-like hypersurface in the usual way is not sufficient to 
uniquely determine its time evolution. In contradistinction to asymptotically flat or de Sitter space-times one also has 
to specify suitable boundary conditions on the time-like conformal boundary of the asymptotically AdS space-time at 
(null and spatial) infinity. 

In Ref. [24 it has been found that the same system, a minimally coupled, free, massless scalar field coupled to gravity, 
as in Ref. |lCll| admits spherically symmetric, time-periodic regular solutions, whose mass is finite. Such localized, 
time-periodic solutions shall be referred to as AdS breathers. In Ref. [2^ AdS breathers have been constructed both 
by perturbative expansion and by direct numerical calculation, although, possibly because of the letter form of the 
paper, concrete results have been presented only for 4-|-1 spacetime dimensions. Both the perturbative and the direct 
numerical methods in are based on the expansion in terms of the eigenfunctions of the linearized problem. In 
case of even spatial dimensions, to each order in the perturbative expansion, all metric and scalar field functions can 
be expressed as a sum of finite number of these terms. This does not hold for odd spatial dimensions, making the 
small-amplitude expansion formalism in not applicable in practice in that case. Similarly, the presence of slowly 
decreasing coefficients in terms of the linearized solutions is likely to make the direct numerical method to determine 
the breather in impractical for even space dimensions. 

Breather-type solutions are of obvious interest and it is a particularly important question to clarify their main 
properties, in particular their stability. In this paper we present a numerical construction of spherically symmetric 
AdS breathers using the Kadath library [2^, [2^ , on the one hand, and compute them analytically by a perturbative 
expansion in the amplitude, on the other hand. Importantly, we also investigate the stability of the AdS breathers 
using a time evolution code. We concentrate on the d = 3 and d = 4 dimensional cases, but our methods work in 
arbitrary dimensions. 

Our numerical construction of the AdS breathers is based on deforming the well known linearized AdS breathers 
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having eigenfrequencies a;„ = d + 2n where n = 0, 1 ,. .. is the node-number, whereby the amplitude of the 
solution becomes frequency dependent. We have mostly constructed breathers corresponding to the nodeless mode, 
n = 0. In Ref. [1^ gravitational geons - localized, asymptotically AdS vacuum solutions of the Einstein equations 
with cosmological constant - have been calculated in 3 -I- 1 dimensions. Nevertheless, to the best of our knowledge, 
scalar field AdS breathers in 3 -I- 1 dimensions have not been constructed numerically in the literature before. In case 
of 4-1-1 dimensional spacetimes, careful comparison of our solutions with those of Ref. [2^ shows very good agreement, 
providing rather nontrivial checks on the validity of both approaches. We note that even for 4-1-1 dimensions our 
approach is able to provide significantly larger amplitude solutions than those presented in . The likely reason for 
this is the choice of the parameter e in [2^ , which reaches a maximal value already before the amplitude of breather 
grows to the value where the mass is maximal. 

Importantly, we have found that AdS breathers whose frequencies lie between Wg and d, where cOg ~ 2.253 for d = 3 
and uJs ~ 3.548 for d = 4, are stable. The ch ang e of stability occurs where the mass of the AdS breather is maximal, 
in analogy to what happens for boson stars |29l - K^ . 

We also hnd a rather complicated resonance structure of the nodeless breathers both in d = 3 and d = 4, where the 
amplitude of certain higher Fourier modes of the scalar field increases abruptly. We have observed several such peaks, 
although here we only present numerical evidence for those which are easiest to find numerically. The detailed study 
of these resonances would require further extremely high precision numerical runs. However, the importance of their 
existence is reduced by the fact that apparently they all lie in the high amplitude unstable domain. The smallest 
frequency where our numerical code converges is around a;~2.1ind = 3 and a;~3.5ind = 4. The existence of a 
minimal frequency for each d and node number, n is expected. 

The paper is organized in the following way. In Sec.|lT]we present the spectral methods used for the construction of 
time-periodic solutions. First, we compute the linearized solutions of the problem, and establish the agreement with 
the analytical results. These solutions are used as initial guess for the iteration procedure for obtaining the solutions 
of the full nonlinear problem in d = 3 and d = 4. For d = 4 spatial dimensions we find very good agreement between 
our results and those of [l^. In Sec. m the numerical time evolution code is presented. We show that numerically 
constructed AdS breathers are stable under time evolution as long as their frequencies are above Wg, and those with 
lower frequencies are unstable, and form black holes. In Sec. IlYla small-amplitude expansion procedure is presented 
for arbitrary dimensions. We have performed the expansion up to fourth order, and excellent agreement with the 
numerical results have been found for as high values of the expansion parameter as e ~ 1. This strongly suggest that 
our perturbative expansion has a finite convergence radius. 


II. CONSTRUCTION OF PERIODIC SOLUTIONS BY SPECTRAL METHOD 
A. Field equations with a massless scalar field 


We consider d -I-1 dimensional Einstein’s equations with a negative cosmological constant A, 

Gfiv + Ag^i^ = SnGTfj^i, , (1) 

where the stress-energy tensor is provided by a minimally coupled zero mass real Klein-Gordon scalar field 




( 2 ) 


where all the Greek indices are d -I- 1-dimensional. The divergence of the Einstein’s equations yields the wave equation 


= 0 . (3) 

Since for d = 1 spatial dimensions both the Einstein tensor and the energy momentum tensor is traceless, equation 
0 can only be satisfied for zero cosmological constant. Hence in the following we assume d > 2. 


B. Choice of coordinates 

The spectral code used in this section works in coordinates such that G = 1. This is to be contrasted with the 
evolution code of Sec. IHII which employed SttG = d — 1. The reader should keep that in mind when comparing results 
coming from various sources. 
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There are many possible choices of coordinates for the AdS spacetime, and it is also the case for spacetimes that 
are only asymptotically AdS. For the spectral methods applied in this section a natural choice is to apply isotropic 
coordinates. In those coordinates the d + 1 dimensional AdS metric reads 


ds2 = - 



( 4 ) 


where p is defined as p = r/L, and dr2^_j^ is the metric on the d—\ dimensional unit sphere. The radius r goes from 
0 to L, so that p G [0,1]. The maximum radius L is related to the cosmological constant by 


d{d- 1) 

2L2 


( 5 ) 


We note that rescaling the time and radial coordinates one can always set L = 1, which we will employ in actual 
numerical calculations. 

In the language of the d + 1 formalism, O means that the spacetime is simply described by a lapse function N, a 
conformal factor and no shift vector, so that 

ds^ = (dr^ + . (6) 

A simple identification yields the expression of the metric potential in the AdS case: 


1 + 




1 — p2 



( 7 ) 


Close to r —>■ L, the two metric fields diverge. More precisely, they behave like 



( 8 ) 


where e = L — r. 

When one deals with spacetimes that are only asymptotically AdS, the expressions ([7]) no longer hold. The metric 
can still be described by the form ([6]). As we are interested in spherically symmetric and periodic solutions, the fields 
N and 'k depend only on r and t. The proper geometry at p = 1 is recovered by imposing the AdS asymptotic 
behaviors ([8]). 

Various d + 1 projections of the stress-energy tensor will appear in the decomposition of the field equations. The 
energy density is given by 




2^4 


and the trace of the spatial part of the stress-energy tensor is 

5 = 




2V2 


(9) 


( 10 ) 


C. d-|-l equations 


Even if there is no shift vector, due to the fact that the fields depend on t, there is an extrinsic curvature tensor. Its 
expression comes from the kinematic definition of the extrinsic curvature tensor, which simply reads: dtjij = —2NKij, 
where 'jij = is the d-dimensional spatial metric and ftj the d-dimensional flat spatial metric. It follows that 


K^J = 


2 

N 



( 11 ) 


The Hamiltonian constraint equation, with a cosmological constant, reads R + — 2A = IGirGE, where 

E is the d -I- 1 energy density given by ([9]) . It leads to 


-h (d - 1) ^-h (d - 3) 


■q, ]V2 i .t) 4^2 


47rG 
d — 1 




( 12 ) 
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The second equation that we solve is obtained by taking the trace of the evolution equation for Kij. For zero shift 
this takes the general form 


= -D^D,N - - 


2N 
d — 1 


A + SnGN 


id-2)E + S 

d^ 


(13) 


and using our variables it gives the following result 

Nrr + {d-l)—+2(d-2)Nr— = 
’ r ’ 4' 

^4 2dN^t-^,t ^ 2d 


N 'i> N N 'i> N \ 'if 


+ N'if'^^ +8TTGN'i/ 


^{d -2)E + S 
d — I 


where S is the trace of the spatial stress-energy tensor given by (nnj. The wave equation ([3]) reads 


(14) 


1 

7^ 


t Nt ' 


1 


i + (d - 1) ^ + ^4>,r + 2 (d - 2) 

r N w 


= 0 


(15) 


D. Equations in the region close to infinity 


We aim at solving the equations numerically, using spectral methods (see Sec. Ill El) . In order to maintain the 
precision of spectral methods, it is then desirable to work with C°° functions. This is not the case of N and 'k that 
diverge at r = L. Motivated by the asymptotics ([51), in the region close to L, the following auxiliary variables are 
used 


n = 


1 

N 


/ = 


1 

i]/2 ■ 


(16) 


Both n and / vanish linearly when r ^ L and are well described by our spectral expansion. Equations (US, (113 
and m can be rewritten in terms of n and / and multiplied by the appropriate powers to get rid of all diverging 
quantities. They give rise (respectively) to the following set of equations: 


, or4/„ ^2 , ,, ^ J- r ^ , 4 / . n2 


-nf n^rr - nf 


-n^r + ‘2-f {n,r) +{d-2)nrn^rf,r + drGff^tt-‘2.d7G{f^t) + 


dn'^fn^tf,t - 


L2 


- SnGn^f^ 


d- 1 


+ dn^f^t4’,t - n^fn^tP,t + nfp^rr + nf' 


.d — 1 


= 0 


P,r - fn^rP,r - (d - 2) nf'^f^rP.r = 0 


(17) 

(18) 
(19) 


The numerical boundary conditions at r = L are given by looking at the terms that vanish the less rapidly in Eqs. 
dm and (jlSp . respectively. It follows that at infinity 


2np + {d-2) n^rf,r = ’ 

fr = -- ■ 

J,r ^ 


( 20 ) 

( 21 ) 


Those conditions are supplemented by the fact that the field vanishes at infinity: (j> {r = L) = 0. It is also possible, 
from the numerical solutions, to verify that the conditions (|2ni) and ([2T1) do lead to asymptotically AdS solutions (i.e. 
solutions for which / = 0 and n = 0 when r = L). 


E. Numerical methods 


The numerical setting used in this paper is very similar to the one employed in our previous work [3^ . The Kadath 
library [2^ [2^ is used to numerically solve the equations at hand. The library enables the use of spectral methods. 
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Here a two-dimensional setting is used, for the fields do depend only on r and t. The radial coordinate ranges from 
0 to L and several numerical domains (typically 4) are used. In each domain, there is an affine law between the true 
radial coordinate r and the numerical one r*. It ensures that r* £ [~1,1] ([0,1] near the origin). Spectral expansion 
is performed with respect to r* and Chebyshev polynomials are used (even ones only near the origin). The influence 
of the number of radial domains on the precision is difficult to assert beforehand. However there is no real necessity 
to measure this effect. What is more meaningful is to verify that the error converges with the number of spectral 
coefficients fast enough that the required accuracy can be achieved in practice. This is the object of Sec. Ill FI 

As far as time is concerned, one aims at finding periodic solutions so that the fields are expanded onto trigonometrical 
functions. More precisely, it can be shown that the scalar field cj) can be expanded onto odd cosines and the metric 
fields onto even ones. Only one temporal domain is used where the numerical time t* = cot ranges from 0 to tt/2. 

The mathematical system to be solved consists of Eqs (HID, (HI and (1151) . In the numerical domain that contains 
r = L they are replaced by Eqs. (ED, (ED and (ED- Regularity of the fields is ensured at r = 0 and the appropriate 
boundary conditions at infinity are discussed in Sec. Ill PI and given by Eqs. (l20l) and (|2T]) and by demanding that 
the field vanishes at r = L. 

The resulting system is solved using a Newton-Raphson iteration. This requires the use of an initial guess configu¬ 
ration. One possible method consists in using a result coming from a small amplitude expansion. Given the fact that 
the stress-energy tensor is quadratic in (j), at lowest order one can fix the geometry to AdS (Eq. (j4D) and solve only for 
the wave equation. In order to avoid the trivial solution </> = 0, we use a similar technique than the one used for the 
study of massive scalar fields on fixed adS background. We impose that the first harmonic takes a given value at an 
intermediate radius r^ud ■ 'k (r'mid) = 0.1*cos (uit). This condition is used as an outer boundary condition for r < Tmid 
and as an inner boundary condition for r > rmid- Doing so one finds a non-zero solution but which radial derivative 
is, in general, not continuous at r = rmid- Only for some discrete values of u is the solution C°°, representing a true 
solution. EigureH] shows some results for the linear solution in the case d = 3. In this case one recovers the fact that 
the (angular) frequencies of the small amplitude solutions are ujn = 3 -I- 2n and that the solution with iOn has n nodes. 
The error on the continuity is of the order of 10“^°, being typically the overall precision of the code (due to finite 
number of coefficients and round-off errors). 




P 

FIG. 1: The first panel shows the error on the radial derivative as a function of uj, for d = 3. The location of the zeros correspond 
to the true solutions of the small amplitude expansion. The second panel shows the different solutions with increasing number 
of nodes, i.e. for increasing values of uj. 

Small amplitude configurations are used as an initial guess for computing non-linear solutions. One works at a 
frequency close to the linear one: w = a;„ — and one scales the linear solution so that ()> (r = 0, t = 0) = 5^. With 
an appropriate choice of and it is usually possible to achieve convergence to a non-vanishing solution of the 
full system. Once a first configuration is found, a sequence can be constructed by slowly varying oj. 


F. Assertion of the precision 

A first check is concerned with the spectral convergence of the solutions. By this one means a convergence faster 
than any power-law, as the resolution increases. In practice an exponential convergence is often observed. This 
is one of the most striking advantage of using spectral methods. One resolution is a set of two integers {Nr,Nt), 














6 


being respectively the number of coefficients used for the radial variable (in each radial domain) and for the time. 
For largely empirical reasons, the following resolutions are considered (keeping roughly Nr ~ ^Nt) : (11,7), (13,9), 
(17,11), (21,13), (25,15) and (33,17). Figure [2] shows the difference between the value of (r = 0) obtained at the 
highest resolution and the one for various lower resolutions (i.e. in a sense considering the value for (33,17) as being 
the ’’true” value). Results are shown for both d = 3 and d = 4. In the case d = 4 the spectral convergence is easily 
seen. For the d = 3 the convergence is actually so fast that the round-off error of 10“^^ is rapidly reached. It is 
true that convergence of the result does not always mean precision but it is usually a rather good indicator. In the 
following, the various results are obtained for the highest resolution, that is for {Nr,Nt) = (33,17). 



FIG. 2: Difference between the result for {Nr, Nt) = (33,17) and various lower resolution results. The y-axis denotes difference 
in the value of the scalar field at the origin and the x-axis the number of radial coefficients Nr- 

In order to further assert the validity of our numerical code, results are compared to the ones obtained, for d = 4 
in [U. It is a strong test because the methods and coordinates are very different. It is indeed so much the case that 
even just comparing the results is not at all straightforward. First one needs to compute, from our solution, the value 
of the parameter e used in [2^. It can be shown to be the scalar product of (j){t = 0) and the lowest order linear 
solution, which is eo in the language of [1^. Using their equation (11), one can see that eo (x) = 2-\/6cos‘* a:. The 
scalar product in question is the one used in [^ . that is the one with a weight function tan^ x. It follows then that 

/.7r/2 

e = 2\/6 / 0 (r, t = 0) cos a; sin^ xdx. (22) 

Jo 

We note that this e is different from the small-amplitude expansion parameter e used in Sec. IIVI of our paper. 

In order to compute the integral appearing in Eq. (l2^ , one needs to express the radial coordinate r as a function of 
the coordinate x used in [2^. A direct comparison of the line elements, for t = 0, shows that r must be the solution 
of 4'^ (r) r = Ltanx, which can be solved easily using, for instance, a secant method. One must also recall that the 
spectral numerical code in this paper uses units such that G = 1 whereas [2^ uses SttG = 3 so that e must be scaled 
by an additional factor y^87r/3. 

The last complication comes from the fact that we measure the frequency ut at infinity whereas the 0 used in [2^ 
is defined at the origin. One can show that the two relates by O = uj/Nq (r = 0) where Nq is the constant part of our 
lapse function (remember that N = Nj cos {2jcjt) ; see Sec. Ill El) . 

Figure [3] shows the comparison between our results and the ones from More precisely, it shows the value of O, 
as a function of e for both codes. The agreement is so good that the results are indistinguishable by eye and a relative 
difference of order 10“^ is measured on the whole range of e. We view this as being sufficiently small to confirm the 
validity both works. More precision (i.e. smaller error) could probably be achieved by increasing the precision of the 
comparison itself (for instance when computing the integral (1221) 1. 

A striking feature of Fig. |3] is the fact that the parameter e has a maximum value of about 0.087. This explains 
why the code used in fails to converge for e > 0.085. Nevertheless, the sequence of solution continues to exist, 
with increasing values of U. This indicates that the variable e is probably not the best choice to parametrize the 
various solutions. In our case, the parameter w is used instead and it seems to lead to an easier exploration of the 
parameter space, as illustrated by the results shown in Sec. Ill Gl In order to facilitate the comparison with the results 
in [ 2 ^ we plot the dependence of the central frequency U on the asymptotic frequency w on Fig. 0] 
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FIG. 3: Values of the central frequency as a function of the expansi on p arameter e used in [^. The solid line denotes the 
results from this paper and the circles the data taken from Table II of . 


d=4 



FIG. 4: Oscillation frequency If measured by a central observer, as function of the frequency u observed by a faraway observer, 
for d = 4 spatial dimensions. For small amplitude states both frequencies approach the value 4. 


G. Numerical results 

In this section, some numerical results are presented for d = 3 and d = 4. We consider only configurations that 
connect to the nodeless linear ones. First, Fig. [5]shows the central value of the first three modes of the scalar field 
(remember that only the odd modes are non-zero). As expected all the various modes go to zero when one approaches 
the frequency of the linear solution (which is w = d in this case). Even in the case d = 4, the curves are relatively 
monotonic and it confirms the fact that the behavior exhibited in Fig. [3] comes solely from the use of the variable e. 

Figure [6] shows the spatial structure of the solution for one particular example, corresponding to d = 3 and w = 2.3. 
The first panel shows the various modes of the scalar field (pn and the second one the values of the metric function N 
and 41, at t = 0. The dashed curves represent the linear solution (first panel) and the metric fields corresponding 
to pure AdS spacetime (second panel). Although with our code the solutions stop at 2.104 for d = 3, and at 3.48 
for d = 4, we do not think that the solutions simply disappear for such values of the frequency, rather we expect 
the families of solutions to continue. It seems plausible, in analogy with boson stars, that close to these points the 
solution curves turn back, and the new branch continues for increasing values of the frequency. 

One of the most interesting findings of our previous work on fields on fixed AdS background was the existence 
of many resonances (see Sec. V-C of [1^). The solutions obtained in this paper are asymptotically AdS and we also 
found a number of resonances in analogy with the fixed AdS background ones. The resonances are, however, relatively 
difficult to find because they are very narrow and noticeable only in rather high order modes. Two such resonances 
are shown in Fig. [71 for d = 3 (first panel) and d = 4 (second panel). In order to check that they are real Fig. [7] also 
shows some results obtained with two lower resolutions (the circles correspond to {Nr^ Nt) = (21,13) and the squares 
to (25,15)). For d = 3 the various resolutions are indistinguishable and a clear convergence is seen in the case d = 4. 




d=3 


d=4 



FIG. 5: Values of the various modes <!>„ at the origin, for d = 3 (first panel) and d = 4 (second panel), as a function of u). The 
modes are scaled for plotting convenience. 



FIG. 6: Solution for d = 3 and w = 2.3. The first panel shows the first modes of the scalar field and the second one the metric 
fields at t = 0. The dashed lines are the linear solution (first panel) and the solution for the exact AdS spacetime (second 
panel). 


This makes us very confident that the observed resonances are genuine and not numerical artifacts. 

We have not seen major differences in the resonances between the cases d = 3 and d = 4, in terms of shape (this 
is to be contrasted with the ones observed in which shape could vary from one case to the other). On Fig. 0 
the ninth mode has been plotted for d = 3 and the fifth one for d=4. They correspond to the modes where the the 
resonances are the easiest to see: they are also visible, although much less noticeable, in the other modes. Let us state 
that in theory the resonance is expected to be present in every mode due to the coupling between them. However, 
this would require a very precise fine-tuning of the frequency to be seen clearly. 

In (3^ a relatively simple explanation for the appearance of the resonances was given: basically, one of the modes 
became dominant and was approaching one of the small amplitude solutions (see the second panel of Fig. 9 in [3^). 
Unfortunately no such simple interpretation is found to hold in the case at hand. It might come from the fact that the 
resonances correspond to cases far from pure AdS background and the metric fields play an important role. Another 
difference with |32| is the fact that the resonances are less numerous, occurring in higher order modes and narrower. 
All those reasons makes them relatively more difficult to find. We are very confident that this is a genuine effect but 
a more detailed study would be required to enhance our comprehension of those resonances. 
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d=3 



d=4 



FIG. 7: Mode (j)n at the origin, around two resonances fonnd for n = 9 and d = 3 (first panel) and n = 5 ad d = 4 (second panel). 
The solid lines denote the high resolntion resnlts {Nr, Nt) = (33,17), the sqnares the middle resolution ones {Nr,Nt) = (25,15), 
and the circles the lower resolution ones {Nr, Nt) = (21,13). 

III. TIME EVOLUTION 
A. Field equations 

In the following sections of the paper we use the same coordinate system and variables as in ©[HI- We look for 
spherically symmetric solutions with metric 

ds^ = -^— ( —Ae~^^di^ + + sin^ a: , (23) 

where L is related to the cosmological constant A by (O, and A and 8 are functions of the coordinates t and x. We 
assume that the center of symmetry is regular at x = 0 and that the metric approaches that of the anti-de Sitter 
one at infinity x = tt/2. With the metric form choice ( 1 ^^ the anti-de Sitter metric is simply A = 1 and 5 = 0. 
The relation with the metric form ([51) follows from = Ltanx, and since in general 'F depends on both t and r, 
the constant r surfaces are different from the constant x ones. Similarly, the constant t and t hypersurfaces are not 
the same either. However, at the moment of time reversal symmetry, the constant t = 0 and t = 0 surfaces agree, 
which makes easier the comparison of the results obtained with the two different coordinate system choices. We note 
however, that the calculation of Fourier components in time, such as those shown on Fig. [B] and [71 is coordinate 
system dependent, and the comparison of Fourier modes in two different systems can be particularly difficult. 

By Birkhoff’s theorem, the only spherically symmetric vacuum solution is the Schwarzschild anti-de Sitter metric, 
which is given by 5 = 0 and 


A = 


1-ms 


cos X 


sin'^ X 


(24) 


where the constant ms is related to the total mass Ms of the Schwarzschild anti-de Sitter spacetime by Ms = jms, 
with 


7 = 


d-1 

IGttG 


Hd-iL 


d-2 


(25) 


and the surface area of the d — 1 dimensional sphere is jT{dj2). 

In general spherically symmetric spacetimes there is a naturally defined radius function R, defined in terms of the 
area of the symmetry spheres. Using this function one can construct the Misner-Sharp energy (or local mass) function, 
which for negative cosmological constant can be defined as 


AI = ym , 


( 26 ) 
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where 


^d-2 




j,2 \ 

I?)- 


(27) 


For the Schwarzschild anti-de Sitter metric M = M 5 , constant. At infinity the mass function M tends to the Abbott- 
Deser mass Mad = I'm ad [l^- In the coordinate system (1231) the radius function is i? = Ltana;, and the function m 
can be expressed by the metric function A as 


sin'^ ^ X 


m = 


cos“ X 


(1-A) . 


(28) 


This shows that in order to have a finite total mass, the function A must tend to 1 at infinity x = 7 r /2 as ( 7 r /2 — x)'^. 
Introducing the variables 


$ = 


dcj) 


n = 


d(j) 


5x ’ "" A di ' 

the wave equation ([3]) in the coordinate system (1231) can be written into the form 

9n 1 d /Atan‘^“^x 


dt tan' 


d-l , 


dx 






The equality of mixed partial derivatives of (j) yields the evolution equation for $, 

9$ d f A^ 
di dx 

The (t, t) and the (x, x) components of the Einstein’s equations give 

85 SttG 


dx 


d- 1 


sin X cos X ($2 +n2) ^ 


dA , 5(5 ., / , 1/1 

= A— —h dtanx + (d — 2) cotx (1 — A) . 
ax ox 


(29) 

(30) 

(31) 

(32) 

(33) 


If $ and n are given on a constant i hypersurface, then the metric variable 5 can be calculated by integrating (1321) . 
and A can be obtained by solving the differential equation (1551) . Equation (15^ shows that it is advantageous to use 
units (or rescale (j)) such that 


SttG = d — 1 , 


(34) 


which we will assume in the next sections, and by this choice we will agree with the notation of m- The (t, x) 
component of the Einstein’s equations give 


dA 

'dt 


IOttG , 

—-sm X cos x—r *1*11 , 

d-l 


(35) 


which is not independent from the previous equations. The remaining angular components of the Einstein’s equations 
turn out to be equivalent to the wave equation (1501) after substituting the derivatives of A and 5 from (1551) . (1551) and 
?51) . 

is not independent of A, but if we define 

(36) 


The spatial derivative of A, as given by the right hand side of 


/i = me , i.e. 


M = 


sm 


d-2 


X 1 — A 


then the expression for its derivative becomes fj, independent [ 2 ^, 

dfi SttG tan*^”^ x 


($2 +n2). 


(37) 


If 5 is already calculated by dSlD, then fj, can be obtained by integrating this equation, and next, A can be calculated 
by a simple algebraic operation using 
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Let us denote the coordinate distance from infinity by y = 7r/2 — cc. It is easy to check, that for a fixed anti-de 
Sitter background, i.e. for A = 1 and <5 = 0, the leading order behavior of a massless scalar field 4> is either or y‘^. It 
can be shown that if a finite mass can be associated with the self gravitating system, i.e. A—1 tends to zero at y = 0 
as y‘^, as in HMD, then the mass is time independent, and (p must necessarily tend to a time independent constant at 
infinity, which can be set to zero. The expansion of the functions then starts as 


P = PdyUOiy‘^+^) , 

(38) 

A = l-mADy‘^ + Oiy‘^+^) , 

(39) 

6 = Soo+0{y^‘^) , 

(40) 


where the constant ruAD is the Abbott-Deser mass of the system, while (pd and doo are functions of t. All higher 
coefficients of the powers of y va (p, A and 5 can be expressed in terms of ttiad, (pd, Soo and their derivatives. For 
even d only even powers of y occur in the expansion of p, A and <5, but for odd d both odd and even powers occur. 
Introducing a new time coordinate as some function of t it is always possible to set 

<5oo = 0 , (41) 

which we assume in the following, as it was also done in [T^, [2^. This means that asymptotically i will give the 
standard anti-de Sitter time coordinate. Alternatively, it would be possible to set d = 0 at the center a: = 0, as it was 
done in [3,113, HI. In that case Li would give the proper time at the center, but the component of the metric 
would be time dependent at large distances from our localized object. In the following we will always assume that 6 
tends to zero at infinity. 

We look for solutions which are regular at the center of symmetry a; = 0, which requires that A tends to I there, 
moreover 6 and p tend to some time dependent finite values. It also follows that in the power series expansions of A, 
(5 and p only even powers of x appear, which implies that the functions have mirror symmetry there. 


B. Time-evolution code 


Our numerical code, which we use for the calculation of the time evolution of spherically symmetric configurations, 
is a fourth order method of lines code very similar to the one employed in , which is described in detail in 

[^ . The coordinate system (1^ is used, and the funtions A, 6, $ and 11 depend on the coordinates i and x. 

We use a uniform grid with grid spacing Ax, and whenever possible we use symmetric fourth-order stencils to 
calculate radial derivatives. This can always be done near the center by taking into account that the functions p^ A 
and 5 are symmetric at x = 0. Near the outer boundary at x = 7r/2 we use the known asymptotic behaviour for the 
calculation of derivatives. If a function goes to zero as / = ?/“(co + C 2 y‘^ + •••), where a is a known positive integer, 
then we calculate cq and C 2 from the value of / at the last two grid points before the boundary, and for the derivative 
we substitute the corresponding value of —f. 

If we need to calculate the radial derivative of a product, such as for example that of Ae“‘^n in the evolution equation 
m for $, then instead of expanding by the product rule of derivatives and calculating separately the derivative of 
n, A and 6, we calculate first the value of the product Ae“‘^n at each grid point and calculate the derivative from 
that. This way it is possible to obtain a more stable numerical method. The equation (1301) can be similarly used in 
order to obtain the time-evolution of 11, but it generates numerical instabilities both at the center and at the outer 
boundary. Hence it is used only for tt/ 8 < x < Stt/S. For x < tt/S and for x > Stt/S we apply the method described 
in [s^ by rewriting (l30l) into the form 


OT_/ d-1 \ d d-1 9 \ 

di V cos(2x) J dx \ J 2 ^ dx\e^ sinx cosx ) 


(42) 


This way we only have to calculate the derivatives of functions which are less singular at the boundaries, and we 
obtain a stable numerical method. In the last term of (HU) the value of $/ sinx at the center x = 0 is substituted 
by the derivative of $ using I’Hopital’s rule. Except for the last two grid points before the boundary at x = 7r/2, 
dissipative terms proportional to the sixth radial derivatives of $ and H are added to the time-evolution equations, 
in order to make sure that the numerical evolution is stable. By choosing both proportionality factors as cr(Ax)®, the 
numerical method remains fourth order, and the value of the constant a can be chosen, for example, to be 0.03. 

Time-evolution of $ and H are calculated at each grid point by a fourth-order Runge-Kutta scheme. At each time 
step, including the Runge-Kutta substeps, 5 is calculated by numerically integrating the right-hand side of (15^ from 
the outer boundary, where <5 = 0. The variable A is calculated by first integrating (l37l) from the center, where y, = 0, 
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and then calculating A using (l36l) . This way we avoid integrating the differential equation (l33l) for A by a Runge-Kutta 
scheme, which would require getting midpoint values at Ax/2 by numerical interpolation. In our method, the forth 
order accurate numerical integration of 6 and is done by fitting third order polynomials at four neighboring grid 
points, and using the integral of this interpolating function between two neighboring grid points. The asymptotic 
behavior of the integrand of S at the outer boundary is / = + C 2 j/^), while the behavior of /r near the center 

is of the form / = + C 2 X^). The constants cq and ci are calculated from the values of the integrand at the 

first two points, and the integral up to these points is represented by the integral of /. The typical number of spatial 
grid points used on our numerical simulations were between 2® and 2^^ points. 


C. Results obtained by the time-evolution code 


We have used several configurations obtained by the spectral numerical code in Sec. |n]as initial data for our time- 
evolution code. We take the field value at the time slice where the configuration is time reversal symmetric, i.e when 
n = 0. Only the field $ is taken from the spectral code, appropriately changing the radial coordinates from r to x. 
Values of A and S are calculated by the same method on the initial time-slice as on all the later time-slices. Firstly, 
this procedure confirms that the time evolution of the initial data indeed leads to time-periodic solutions with the 
expected oscillation frequency, precise to several digits. Secondly, it enables us to study the stability properties of 
these configurations. 

Small amplitude configurations turn out to be always stable. These have frequency w, which is slightly lower than 
the frequency ojq = d of the nodeless linearized solution. These remain oscillating with essentially constant frequency in 
the time evolution code for arbitrary long time-periods. For higher amplitude configurations, the oscillation frequency 
is lower, but the behavior remains very similar. In order to illustrate this, on Fig. we show the time evolution of the 
frequency when numerically evolving initial data corresponding to w = 2.3, in case of d = 3, with increasing resolutions 
in the time-evolution code. The frequency is obtained by calculating the time difference between the current and the 




FIG. 8: Difference from the expected frequency of the oscillations as a function of time, evolving initial data corresponding to 
uj = 2.3 and d = 3. An initial and a much later period is shown for numerical resolutions from 2^ to 2^^ points . 


previous maximums of the oscillating scalar field at the center. 

However, as the amplitude increases and the oscillation frequency becomes lower, below a certain frequency Wg, 
the configurations become unstable. For d = 3 spatial dimensions oJs = 2.253, and for d = 4 the instability starts 
below ujg = 3.548. In both cases this corresponds to the frequency where the mass of the configuration is maximal, 
as it is usually the case for physically relevant localized objects. We note however, that for self interacting scalar 
fields on a fixed AdS background instability generally occurs already before the energy of the configuration reaches 
a maximum [^. On Fig. we show the frequency dependence of the mass ruAD in the d = 4 dimensional case. The 
mass is calculated from the results of the spectral code using (l?7l) at some large radius. A similar figure for d = 3 will 
be given in Sec. IIV Gl together with the corresponding small-amplitude expansion results. 

In the d = 4 case Ws = 2.253 corresponds to the frequency Og = 6.325 as measured by a central observer (see 
Fig. 0]). The configurations represented by the circles on Fig. [31 which were taken from [^, are all in the stable 
domain, instability only occurs for higher amplitude configurations. We note that the instability does not appear at 
the turning point where e is maximal on Fig. [31 it begins at somewhat higher H. 
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FIG. 9: Frequency dependence of the mass of the configurations for d = 4 spatial dimensions. 


The narrow resonances that we have found using the spectral numerical code in Sec. Ill Gl are all in the unstable 
domain, both for d = 3 and d = 4. 

The unstable solutions initially oscillate with the expected frequency, however after about a few dozen oscillations 
the oscillation amplitude starts changing, and eventually the scalar field collapses into a black hole. The metric 
function A can be seen to approach zero at some radius x, which signals the appearance of an event horizon. The 
more precisely the initial data is given, the longer the configuration oscillates, and the later the collapse occurs, which 
indicates that very likely there is a perturbation mode which makes these high-amplitude periodic solutions unstable. 

As a concrete example, on Fig. [TUlwe show the time evolution of the energy density 


^ ACOS^X ,n o. 


(43) 


at two places for initial data belonging to the frequency uj = 2.104. The upper plot shows E at the 
center, while the lower one at a radius which is outside the sphere where the event horizon first appears. 


symmetry 
For about 




FIG. 10: Time evolution of the energy density E for an initial data which is in the large amplitude unstable domain, corre¬ 
sponding to ca = 2.104. The upper panel shows E at the center, while the lower one is at a radius that is larger than the one 
where the horizon forms. 

six oscillations the evolution appears to be periodic, but then the unstable mode grows large enough to destroy the 
periodicity, and the scalar field quickly collapses into a black hole. The lower panel shows that the energy density tends 
to zero outside the horizon, while the upper panel shows that it approaches a constant value at the center. However, 
this does not mean that the energy density actually stops increasing at the center, it is only the time foliation which 
slows down very much, because the metric function 5 increases to large values in the central region. Our coordinate 
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system, which is based on a constant area-radius foliation, breaks down when an apparent horizon forms. The specific 
way of how a conhguration collapses to a black hole, and after how many oscillations, depends strongly on what kind 
of small perturbations are already present in the initial data, and also on how strong is the instability for the given 
frequency. The large amplitude configuration which we show on Fig. [TO] is relatively short living even if the initial 
data is calculated with very high precision by the spectral code. 


IV. SMALL-AMPLITUDE EXPANSION 
A. Linearized equations 

We expand around the vacuum anti-de Sitter metric using the coordinate system ([23|) . We use our previously 
developed perturbative framework for oscillons [H, [s^ and oscillatons [13, . In terms of a small parameter £ we 

write 


OO OO OO 

n=l n—1 n—1 


(44) 


where and 5*-"^ are functions of i and x. We substitute this expansion into the wave equation (1301) . and into 

the Einstein’s equations components ( 15 ^ and ( 1551 ) . Since in many cases it is easier to solve, we also use ( 1551 ) . but we 
cannot eliminate the use of (1331) completely. It will turn out later, that the expansion of </> contains only odd powers 
of e, while the expansion of A and S only even powers of s. Assuming that there is only a single one-parameter family 
of solutions that we consider, this follows from the (j) —>■ —(^ symmetry of the system. We also note that our expansion 
parameter e is not the same as the parameter e used in [24| . although to leading order they are proportional to each 
other. 

Since the Einstein’s equations contain (j) quadratically, to first order in e the equations for A and <5 are the same as 
in vacuum. Since the only spherically symmetric vacuum solution is the Schwarzschild anti-de Sitter metric, which 
only has a regular center if the mass is zero, it follows that A^^'^ = 0 and <5^^^ = 0. From the wave equation (|30|1 it 
follows that the leading order coefficient of the scalar field satisfies the differential equation 

dP dx"^ sin x cos x dx 

Let us look for periodically oscillating solutions in the form = pcos{uji), where p only depends on x. Since we 
assume that 6 tends to zero at inhnity, the time coordinate t is the standard anti-de Sitter time, and the frequency 
oj represents the frequency measured by some faraway static observer. However, since the proper time at the center 
of the vacuum anti-de Sitter spacetime in the coordinate system (1231) is not t but Lt, the physical frequency of the 
oscillating solutions is actually not uj but oj/L. Then p has to satisfy the ordinary differential equation 


d^p 

dx^ 


d-1 


sin X cos X dx 




(46) 


Solutions, which are regular at the center and localized in the sense that p tends to zero at infinity, only exist for 
certain discrete values of the frequency, uj = ujm where 


ujn = d + 2n , 

for n > 0 integers. The corresponding regular localized solutions can be given in the form 


Pn = 


{d/2), 


■ cos'^ X 1.V2) (cos(2x)) 


(47) 


(48) 


where P denotes the Jacobi polynomial, and {a)n = a{a -I- I)... (a -I- n — I) denotes the Pochhammer symbol, with 
(q;)o = I. The solutions p„ are normalized in order to make their central values I at x = 0. The integer n labels the 
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number of nodes of the solutions. The first few functions are 


cos'^x 


Po = cos'^ X 

Pi = 

P2 = 

P3 = 


d 


cos'^ X 

[{d 

2d 

cos"^ 

X 


M{d + 2) 


[(d + 1) cos(2a;) — 1], 

[{d + 3) cos(4a;) — 4cos(2a;) + d + 1], 

{(d + 3) [(d + 5) cos(6x) — 6cos(4a;) + 3(d + 1) cos(2a:)] — 6(d + 1)} . 


(49) 

(50) 

(51) 

(52) 


We note that for odd d it may be possible to write (l48l) in terms of trigonometric functions without using orthogonal 
polynomials. For example, for d = 3 


1 


Pn = 


4 sin a; 


sin[2(l + n)x] sin[2(2 + n)x 


I + n 


+ 


2 + n 


(53) 


Rescaling the polynomials as 


2yy(^n d — 1)! 


(54) 


the functions Cn are normalized in order to make them an orthonormal base on the Hilbert space T^([0,7r/2], tan^“^ dx). 

To the linear approximation, the general localized regular solution can be written as a sum of these solutions with 
arbitrary amplitudes Un and phases bn, 


(j)='^an COs{ujJ + bn)Pn ■ 
n—0 


(55) 


For the linearized problem all these solutions give stable configurations, and since all frequencies are integers, the 
whole solution is periodic with frequency a; = 1 or 2. To the linear approximation, the frequency of each mode 
is independent of its amplitude. However, nonlinear effects will make the frequencies amplitude dependent. In the 
small-amplitude expansion procedure this means that each ujn becomes e dependent. In general, one expects that each 
LOn will change independently, becoming unrelated real numbers, and therefore the solution will stop being periodic. 
In Refs. 0 and (l^ it was argued that proceeding to higher orders in the e expansion secularly growing terms, of 
the type tcos(wt), appear in all cases when two or more, but only a finite number of a„ are nonzero. According to 
[l^ in the case when an infinite number of a„ are different from zero, it may be possible to cancel all secular terms. 
In that case however, the solution will be quasiperiodic, consisting of infinite number of modes with not rationally 
related frequencies [^ . 

Since periodic solutions are of a great interest in their own, in this paper we concentrate on constructing those. If 
only a single a„ is nonzero, then all secular terms can be absorbed by changing a;„ as a function of e. This means that 
there is a one-parameter family of nonlinear solutions emerging from each pn mode, and their frequency gradually 
becomes more and more different from as the amplitude grows. We will be mainly interested in the family of 
solutions starting from the nodeless mode po, since that one is expected to be the most stable, and also the one with 
the lowest mass. 


B. Periodic solutions 


In this paper we concentrate on periodic solutions, having a definite frequency ui, which generally depends on the 
oscillation amplitude. In order to include this periodicity in our formalism more conveniently, we introduce a new 
time coordinate 


T = cot . (56) 

In terms of r the oscillation frequencies of all Fourier modes will be integers. Using r each time derivative in the field 
equations (1501) - (I55|) will bring an oj factor, and there will be terms in the equations. We take into account the e 
dependence of w by expanding it as 


U) 


= -<">h + E 

i=i 




(57) 
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where and are some numbers. It will turn out later, that = 0 for odd j. Assuming that there is only a 
one-parameter family of solutions that we consider, this also follows from the (j) —>■ —(f) symmetry of the system. 

In the zero amplitude limit u> should tend to one of the a;„ frequencies, hence = lOn for some n. From now on 
we concentrate on the solution emerging from the nodeless mode, so we assume that 

= ujQ = d . (58) 

We Fourier decompose the coefficients in the small-amplitude expansion of cj), A and S as 


cos (jr) 
j=o 


j=0 


(n) 


cos(jT) 


i=o 


cos(jT) 


(59) 


where ^cid are functions depending on x. We could also have included sin(jT) terms in the expansion, 

but because of the time reversal symmetry of the problem, those terms would turn out to be zero anyway. It will also 
turn out later that = 0 for even j, while = 0 and = 0 for odd j. This is closely related to the 
symmetry. 

Substituting the expansion into the wave equation dani), the equations for all take the form 


dx2 


d-1 

sin X cos X dx 




'2 J,(") 




(60) 


where depends only on the lower e order functions, with I < n. For the first Fourier modes, i.e. for 

j = 1 , the homogeneous part of the equation is always solved by = cos'^x, consequently one can always add a 
constant times cos'^ x to a given inhomogeneous solution. We use this freedom to specify how we parametrize the 
various states by the parameter e. We set 

= latx = 0 , t = 0 , (61) 

= 0 atx = 0 , t = 0 forn> 2 . (62) 


By this choice the parameter e gives the central value of the scalar field 4> at the time when the configuration is time 
reversal symmetric, i.e. at t = 0. We note that our parameter e is different from the parameter e used in [^ . which 
corresponds to the scalar product of ^ at t = 0 with the nodeless linear solution cq. Our e is a monotonous function 
of the frequency in all the domain where we could calculate periodic states with the spectral numerical code, both in 
d = 3 and d = 4 dimensions. 


C. First order in the small-amplitude expansion 

To the leading order in the small-amplitude expansion, i.e. for n = 1, obviously all in (ISOl) are zero. Comparing 
to (I46L we can see that the regular localized solution for is an arbitrary constant times pq. However, other 
may also have localized regular solutions after setting the basic frequency 0 ;^°^ = d. According to (ITTl) . this 
happens if the frequency jd is equal to d -|- 2n for some n > 0 integer, i.e. if (j — l)d = 2n. If the number of spatial 
dimensions d is even, then this holds for all j > 1, and for odd d there are solutions when j is an odd number. If they 
exist, the regular localized solutions for are given by multiples of pn, and have (j — l)d/2 nodes. The sum of all 
these solutions, with arbitrary phase shifts and amplitudes, are valid solutions of the linearized problem. However, 
continuing to higher order in the e expansion, assuming that is nonzero, it turns out that at order there are 
localized regular solutions only if 

=0 if J ^ 1 . (63) 

We have checked this for d = 3, but we expect that it holds for any dimensions. For example, if we allow that 
= cipi COST -I- C 3 P 3 cos(3t) for some constants ci and C 3 in the d = 3 case, then at order it turns out that for 

there are localized regular solutions only if C 1 C 3 =0. The third order secular terms have been explicitly calculated 
in Ref. 0, where it was noted that most of the potentially resonant terms do not in fact arise. In the language of 
[^, the solution is quasiperiodic if the renormalized amplitudes remain constant. Periodicity with a single frequency 
requires further restrictions on the renormalized phases and the resulting frequencies. It is a promising strategy to 
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find a proof, using the methods of Ref. 0 , to show that exact periodicity is possible if and only if there is one single 
mode present in the linear order of the expansion, for arbitrary dimensions. This problem is beyond the scope of the 
present paper and is left for future research. From now on we assume that (16311 holds. Since we plan to interpret e as 
the central amplitude of ^ at t = 0, according to we take 

= cos'^ X . (64) 

Since to first order the metric components and the scalar held decouple, the only regular solution is the anti-de 
Sitter metric, so necessarily 

= 0 , = 0 . (65) 

Proceeding to higher orders in the small-amplitude expansion, it will turn out similarly, that 

= 0 , = 0 for odd n . (66) 


D. Second order in the small-amplitude expansion: metric variables 


Equation ((5^ gives the derivatives of ' in terms of lower order quantities, which can generally be integrated easily. 
Since we would like t to give the anti-de Sitter coordinate time for faraway static observers, we set the integration 
constants in order to make = 0 at a; = 7r/2. The only nonvanishing components at order are 


^ cos^^' X , 


t cos^'^a; [1 — dcos{ 2 x)]. 

^ 4{d+l) ' 

We remind the reader that we use the choice SttG = d — 1 in this section of the paper. 

Equation (I55|) gives algebraic equations for all Aj , except that it does not give any relation for Aq 
nonvanishing component for j > 0 is 


(67) 

( 68 ) 


The only 


^ 2 ^^ = - sin^ x cos^‘^ X 
^ 2 


(69) 


(o's (2'\ ( 2^ 

For even d, the functions Sq , 62 and can be written as sums of hnite number of basis functions. However, 
for odd d inhnite number of Pn is needed, and the convergence is very slow. The reason for this can be seen by power 
series expanding at infinity x = -k 12. The expansion of each contains only {x — 7r/2)‘^“''^^ terms, where k > 0 
integer. For odd d these are all odd powers. However, the expansion of and contain (x — 7r/2)^'^+^^ 

terms which are all even powers, even for odd d. 

After setting the components of A^^^ according to (I69p and setting = 0 for all other j > 0, the field equation 
(l33l) only yields one first order linear differential equation. 


dAp^^ ^ d—2 cos^x^(2) 

J I ^ ^0 

ax sm X cos x 


LL . _ 1 

— sm X cos a; = 0 . 
2 


This can be rewritten into the form 


_d_ 

dx 


sin'^ ^x 


cos“ X 




d-l 


( 2 x) = 0 , 


(70) 


(71) 


so the problem is reduced to integrating integer powers of the sine function. For general d the integral can be given 
in terms of hypergeometric functions, yielding 






COS X 


2d+i 


sin 


cos(2a:) 2 F 1 ( 1 - cos^(2a:) 


2r(^) 


(72) 


( 2 ) 

The constant of integration was set in order to make Ap regular at the center a; = 0. If the number of spatial 
dimensions d is even, then the second argument of the hypergeometric function becomes a negative integer, and the 
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hypergeometric series stops at finite order. In that case the hypergeometric function can be written in terms of a 
Gegenbauer polynomial. For even d 


A 


( 2 ) _ 
0 “ 


(I) COS'^ X 

2d+2r sin'^-^a; 


1 + (cos(2a;)) 


(73) 


The first few specific values for given even dimensions are 


d = 2 : = — cos^ ccsin^ x , (74) 

2 

d = 4: = — - cos"'xsin^ [2 + cos(2a;)] , (75) 

O 

d = 6: = — — cos® a;sin^ x [19 + 18cos(2a;) + 3 cos(4x)] . (76) 

For odd d it is possible to use identities to convert the hypergeometric function in (1721) into a sum of finite number 
of terms, but the result will not be a polynomial in since and cosx. Instead of that approach, it is easier to apply 
directly the expression 2.513 1 in [s^ for the integral of even powers of the sine function, giving for odd d 


^ (2) cos"^ X 


2 / 


fc =0 


d — 1\ sin[2(d — 1 — 2k)x] 
k ) d — 1 — 2k 


The first few specific values for odd spatial dimensions are 

. „ .(2) 9cos®cc., . ,, ,, 

d = 3 : Alf’ = -[4x — sm(4cc) , 

° 64sincr' ^ 

d = 5 : A^'^ = — X ^ 8sin(4x) + sin(8x)] , 

° 2048sin®a:' ' 

. . zlO 'T* 

d = 7 : A}?' = -=—[120cc — 45sin(4cc) + 9sin(8a;) — sin(12a:)] . 

° 49152 sin® CE 


(77) 


(78) 

(79) 

(80) 


(n \ 

For even d the metric function Aq ' can again be expressed as a finite sum of Pn basis functions, but for odd d 

(2) 

infinite number of Pn are required. In this case the power series expansion of Aq for odd d contains both odd and 
even powers of x — 7r/2 at infinity. 


E. Second order in the small-amplitnde expansion: scalar field 


To order the scalar field equation (1501) gives the same form of equations for as (1001) . and the only nonvanishing 
source term is 

=-2d2a;(')cos‘^cc . (81) 

(n) . (n) 

From now on we take (jr^ = G for all those n and j for which F) "^ = 0. Even if there are regular localized solutions 
of the homogeneous problem, we take the trivial zero solution if it is possible. If we would not do this, we expect that 
there would be no localized and regular solutions for some quantities at two orders higher in the e expansion. We 
cannot give a general proof for this, but we have checked it for some concrete n, j and d values. 

In order to calculate the solution of the inhomogeneous equation (l60l) . we need two linearly independent solutions 
of the homogeneous problem. One homogeneous solution for is 

Ml = cos‘^X2Fi ^^(1 + j), ^(1 - j); 1 + , (82) 


which is always localized in the sense that it goes to zero at infinity as {tt/ 2 — x)‘^. When j > 1 and d(j — I) is even, 
the hypergeometric series closes at finite order, and mi becomes proportional to one of the localized regular solutions 
defined in (H51) 


(- 1 )^ 


-Pn 


Ml = 


J 


n = ^(j - 1) ■ 


(83) 
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When d is even, this happens for all j > 1, while for odd d there are localized solutions if and only if j is odd. This 
is important, since it turns out that at each order in e there are nonzero source terms in (1601) only for odd j. 
For odd d a second solution can be given as 


U2 = 2F1 


dj dj d 2 I 

2 ’ 2 ’ 2 ’ J 


sin'^ X 


iFi 


1 - ^(j + 1),1 + ^(j - 1);1 - ^icos^a; 


(84) 


This function is not defined when the third argument of the hypergeometric function, 1 — d/2, is a nonpositive integer. 
When both d and j are odd, the second form becomes a polynomial, 

W 2 =-7- — rf _2 (cos(2a:)) , n=|(j + l)-l, (85) 

where P denotes the Jacobi polynomial and ()„ the Pochhammer symbol. For concrete dimensions it may be possible 
to find simpler expressions. For example for d = 3 and for odd j 

U2 = (—[3j cos{3jx) cot X + sin(3ja;)] . (86) 


For d = 5 and for odd j 


U 2 = 


(_l)(i+i)/2 


5j 


sin^ X 


cos(3a;) cos(5jx) + (25j^ cot^ x — 3) sin(5jx) 


(87) 


The general solution of the inhomogeneous differential equation (1601) can be given in terms of two independent 
solutions of the homogeneous problem, ui and U 2 , as 


i(n) f ^1 T7i(«) J f ^2 T7i(n) 1 


where xi and X 2 are some constants, and the Wronskian is 

du2 dui 

W = Ml—-M2-|— . 

dx dx 

For odd d, when M 2 can be defined by (IMl) . the Wronskian can be calculated to be 

W = dcot*^”^ X . 


( 88 ) 


(89) 


(90) 


By Abel’s differential equation identity, the Wronskian of any two solutions of the homogeneous equation is a constant 
times cot'^”^ x, and hence by appropriately scaling the second solution M 2 we can make (|90l) valid for even d as well. 


For example, for d = 2 we can choose the second solution as 

= 1 — 2 cos^ X log(cot x) , (91) 

1 3 

U 2 ~^'^ = - [7 + 18 cos(2x) + 15 cos(4x)] — - cos^ x [3 — 4 cos(2x) + 5 cos(4x)] log(cot x) , (92) 

and for d = 4 we can choose 

2 

U 2 ~^^ =6 + 3 cos(2x)- 2 -^ log(cot x) , (93) 

sin X 

1 T 

M+^^ = —2 - ^ cos(2x) + 75 cos(6x) — 360 cos(8x) + 495 cos(lOx)] (94) 

sin X 128 

105 

—— cos"^ X [35 — 40 cos(2x) + 60 cos(4x) — 24 cos(6x) + 33 cos(8x)] log(cot x) . (95) 


However, for even d there is a simpler method to solve equation dOHll . without the calculation of the integrals in (1881) . 
In that case the solution can be sought of as a sum of finite numbers of functions, as it was done in [2^. While 
the method in works only in case of even spatial dimensions, our expansion can be performed for arbitrary d. 

If the first solution ui is regular and localized then the other solution U 2 is necessarily singular at the center x = 0 
and non-localized at infinity x = tt/2. In this case, assuming that the source term Fj^ behaves well in (ISOl) . the 
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second term in the inhomogeneous solution (l88l) is always regular and localized, but the first term in (l88l) can be 
regular and localized only if the integral in it vanishes both at a; = 0 and x = 7r/2. The choice Xi = 0 makes the 
integral zero at the center, but in order to make it vanishing at infinity we have to require that 



^ p(") 

W ^ 


dx = 0 . 


(96) 


This gives a condition on Fj"'^ in each case when ui is regular and localized. For example, for j = 1 at order 
is given by dni), and the condition (l%t yields 


from which follows that = 0. For each positive integer n, at order in the small-amplitude expansion the same 
source terms arise for j = 1 as in (jSTl) . with replaced by from which follows that 


= 0 for odd k . 


(98) 


After this, at order all source terms are zero, and we take the trivial solution = 0 for all Fourier 

components. Then = 0, and as we have mentioned earlier, only odd powers of £ remain in the expansion (1441) 
of 4). 


F. Third order in the small-amplitude expansion 


To order there are no scalar field source terms in the equations determining A and S, hence only the trivial 
(S') (3) 

solution is regular, = 0. The same holds to all odd orders in e. The only nonzero source terms in the 

scalar field equations (l60t at order are 


F® =d‘^cos‘^x 
F^^^ = cos^‘^-^ X 


\ a cos^ X J 

(P P /n 7 2 

-h —(3d-h2)cos^x 
4 4 


2cc(2) - 

P{3d 


2d 


2(d -|- 1) 



cos'^ X 


1 -I- cos(2x) 


COS^ X J 


(99) 

( 100 ) 


The integral condition (IMl) automatically holds for F^^\ but for F^^'^ the condition determines the value of The 
integral of the term containing Aq ' can be calculated by integration by parts using (inD, resulting in 


(2) _ (3(i -I- l)\4n F -I- l) 

22(rf+i)r(f)F(d-H|) 


( 101 ) 


In Table U the value of is listed for several values of d. 


The values of and 


can be calculated relatively 


d 

2 

3 

4 

5 

6 

7 

8 

9 

p^^ 

7 

45 

13 

65 

57 

24871 

25 

21735 

20 

128 

42 

256 

286 

16,3840 

221 

262144 


TABLE I: The value of for various spatial dimensions d. 


easily by the integral formula (j88ll for any concrete choice of d, but we could not find a simple expression which is 
valid for arbitrary dimensions. 


G. Higher order results for d = 3 


For d = 3 spatial dimensions we obtain 

,f 3 'i 27 9 Q 27 cos'^x , 

cos^ X -I- —X— —— (cos(2x) — 3) 


-h 


32 
3 cos^ X 
2895872 


64 sin X 

(743109 — 43632 cos(2x) — 1212cos(4x) — 6464cos(6x) -I- 303cos(8x)), 


( 102 ) 
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and 



3 cos^ X 
413696 


(6845 + 4008 cos(2x) + 7692 cos(4a:) 


1368cos(6a:) + 303cos(8a:)). 


(103) 


O') O') 

The solutions of the differential equations determining (p\ and ^3 ' are not unique, since to any solution we can add 
a constant times the solution of the homogeneous equation, which is regular and localized in this case. Namely, one 
can replace by + cipi, and by + C 3 P 3 , where ci and C 3 are some constants, and Pn are the functions 
defined in (l48l) . The constant ci can be expressed by C 3 if we require that at x = 0 we have = 0, which 

ensures that e really gives the central amplitude of 0 at f = 0 to order in the expansion. The constant C 3 can still 
be arbitrary if we perform the expansion only up to order. However, at order it turns out that the equation 
for can have a regular localized solution only if the constant C 3 takes a certain value. The expressions (I102|) and 
(I103|) already give the corrected functions, for which ci = C 3 = 0. However, we note that it is not possible to guess 
the proper functions just by solving the order equations. The solution (I103p does not correspond to a trivial choice 
of the integration constant X 2 in (1551) . such as 0 or tt/2 . 

At £° order, from the condition of the existence of an appropriate solution for (pi , we get the fourth order correction 
to the frequency, which for d = 3 takes the value 


( 4 ) _ 405 _2 1208433249 

“ 1024^ ~ 324337664 


0.177656 . 


(104) 


On Fig. [TT] we compare the small-amplitude expansion results to the precise results calculated by the spectral 
numerical code, by showing the frequency to as the function of e. The expansion gives reasonably good results in 
almost the whole range where the configurations are stable according to the time-evolution code. With the convention 
8 ttG = d — 1 used in this section, e is equivalent to the central value of (p at the moment of time reversal symmetry. 
However, when comparing to the results in Sec. m one should rescale the scalar field with 2 y/Tr, since the spectral 
numerical code was implemented with the choice G = 1. 



FIG. 11: The oscillation frequency ui is shown as a function of the expansion parameter e. For comparison, the small-amplitude 
expansion results, valid to and e"* order are also given. According to the time-evolution code, the configurations with 
frequency below cOs = 2.253 are unstable, so they are represented by a dotted line in that region. 


Another interesting physical quantity is the frequency H as observed by a central observer. This can be calculated 
by integrating the proper time at the center during one oscillation period. From the numerical results of the spectral 
code this can be obtained as H = a;/Ao(r = 0). Using the small-amplitude formalism one has to use the central values 
of and , and the expansion in e results in 

n = (1 -b q(4)^4 _ (105) 


Here = + d/4, which for d = 3 spatial dimensions give = 51/128. For d = 3 




405 2 

-TT^ 

8192 


129387681 


0.398438 . 


(106) 


324337664 
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e 

FIG. 12: Central frequency as a function of the parameter e. The solid line shows our numerical results, while the dashed 
lines give the small-amplitude expansion results. 


On Fig.[T2]we compare the precise numerical value of O with the results obtained from the small-amplitude expansion. 
The agreement with the order result is remarkably good in all the domain for which we could compute configurations 
with the spectral numerical code. The likely reason for this good agreement is that both e and O are taken at the 
center of symmetry. 

The mass of the configurations can be calculated using the expansion (1311). Only the time independent terms 
give contribution to the mass, and the result can be expanded as 

mAD = rn2e‘^ + rme'^ -I- ... . (107) 


For all dimensions, m 2 = 97r/32, and for d = 3 we get 


81 


638469 


rrii = 


TT a; -1.05316 . 


(108) 


2048 11583488 

On Fig.[T3]we compare the precise numerical results to the first two nontrivial orders of the small-amplitude expansion. 



FIG. 13: Mass of the configurations as a function of their frequency for d = 3 spatial dimensions. The and order results 
of the small-amplitude expansion are also shown for comparison. 


It is quite remarkable that the perturbative expansion in the amplitude show excellent agreement with the numerical 
results up to e ~ 1/2. We recall that e determines the value of the amplitude of the field at the center at t = 0, 
therefore it is a physically meaningful quantity. Based on this, we conjecture that the perturbative expansion has a 
large (0(1)) radius of convergence. 
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V. CONCLUSIONS 

In this paper we have constructed spherically symmetric breathers in a massless scalar field theory coupled to 
Einstein’s gravity in d spatial dimensions imposing anti de Sitter asymptotics on space-time. Using a time evolution 
code we have established stability for families of solutions of AdS breathers. We have worked out a perturbative 
constuction of the AdS breathers in the small-amplitude limit for any d, and we conjecture that it has a finite radius 
of convergence. 
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